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ABSTRACT 

We present ground-based observations of the transiting Neptune-mass planet Gl 436b obtained with 
the 3.5-meter telescope at Apache Point Observatory and other supporting telescopes. Included in this 
is an observed transit in early 2005, over two years before the earliest reported transit detection. We 
have compiled all available transit data to date and perform a uniform modeling using the JKTEBOP 
code. We do not detect any transit timing variations of amplitude greater than ^1 minute over the 
~3.3 year baseline. We do however find possible evidence for a self-consistent trend of increasing 
orbital inclination, transit width, and transit depth, which supports the supposition that Gl 436b is 
being perturbed by another planet of < 12 M® in a non-resonant orbit. 
Subject headings: planetary systems — stars: individual (Gliese 436) 



1. INTRODUCTION 

GHese 436 is an M-dwarf (M2.5V) with a mass of 0.45 
M0 and hosts the extrasolar planet Gl 436b, which is 
currently the least massive transiting planet with a mass 
of 23.17 Me ()Torresll200"7[ l. and the only planet known 
to transit an M dwarf. Gl 436b was firs t discov e red vi a 
radial- velocity (RV) variations by Butl er et al.l (|2004[ ). 
who also searched for a photometric transit, but failed 
to detect any si gnal greater than .4%. It was thus 
a surprise when iGillon et all (|2007b[ ) reported the de- 
tection of a transit with a d epth of 0.7% , implying a 
planetary radius of 4.22 R^ (iTorresI [20071 ) and thus a 
compositi on similar to Uranu s and Neptune. In addi - 
tion, both iPeming et al.l ()2007D and lManess et al.l ()2007D 
calculated that the significant eccentricity of the orbit, 
e = 0.15, coupled with its short period of ^2.6 days, 
should result in circularization timescales of ^lO'* years, 
which contrasts with the old age of the system at >6x 10^ 
years. The existence of one or more additional planets 
in the system could be responsible for perturbations to 
Gl 436b's orbit, and thus result in the observed pecu- 
liarities. We consid ered this possibility right after the 
initial publication of lGillon et al.l (|2007b[ ) , and began an 
intensive campaign to observe the photometric transits 
of Gl 436b in order to search for variations in dicative of 
orbital perturbations ('Stringfellow e t a"ni2008[ ). 

Early this year,,Ribas ct al. (2008a1) reported the pos- 
sible detection of a ~5 Me companion in the Gl 436 
system located near the outer 2:1 resonance of Gl 436b 
via analysis of all the RV data compiled to date. Theo- 
retically this planet would be perturbing Gl 436b so as to 
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increase its orbital inclination at a rate of '^0.1 deg yr~^, 
and thus its transit depth and length, so that the non- 
detect ion by 'Butler et alj ()2004f ) and the observed transit 
of G illon e t al. (2007b) were compatible. Since the RV 
detection of this second plan et had a sig n ificant false- 
alarm probability of ~20%, IRibas et all ()2008aD pro- 
posed that confirmation could be achieved through 2008 
observations of Gl 436b's transits, which would show a 
lengthening of trans it dura tion by ~ 2 minutes compared 
to the Gil lon et al I (j2007bl ) data. As well, transit-timing 
variations (TTVs) of several minutes should also be de- 
tectable by observing a sig n ifican t number of transits. 

Recently, lAlonso et al.l (|2008D reported a lack of 
observed inclination changes and TTV evidence for 
the second planet, based on a comparison of a single 
H band light curve obtained in March 2 008 to SAtni 
data taken with Spitzer 2 54 days earlier (jGillon et al.l 
I2007at iDeming et al]|2007[ ). This result, combined with 
additional ra dial velocity measurements (jHoward 1 120081 : 
iBonfild [200 81 that contr adicted the p roposed period of 
the second planet, drove [Ribas et aiT (.2008b. ) to retract 
their claim of the com panion at lAU Symp osium 253. 
However, very recently IShporer" et al.l (|2008D presented 
multiple light curves obtained in May 2007, and could 
not rule out TTVs on the order of a minute. While 
the planet specifically pro posed bv Ribas et al.l (|2008aD 
most likely does not exist, IRibas et al.l ( 2008bD makes a 
strong case that a second planet is still needed to explain 
the peculiarities of Gl 436b, and most likely exists in 
a non-resonant configuration where no strong TTVs 
are induced. Amateur astronomers have been diligent 
in observing Gl 436b since it's initial transit discovery, 
and thus along with this data, published data, and our 
own data, we are able to present a thorough analysis 
of the TTVs, inclination, duration, and depth of the 
transit changes in the Gliese 436 system. We present 
our observations in §2, our modeling and derivation of 
parameters in §3, and explore the observed TTVs and 
parameters of the system over time in §4. 
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Fig. 1. — Top: The V band light curves observed with the APO 
3.5-meter with model fits. Middle: The 3.5-meter data combined, 
phased, and binned in increments of 0.0005 phase. Bottom: The 
transit observed by the NMSU 1-meter telescope on the night of 
January 11 2005 UT. A 3-sigma clip has been applied, and is shown 
with a model fit for which the radii were fixed. 

2. OBSERVATIONS 

We observed Gl 436 (llh42mlls, -f26°42'24" J2000) 
in the V filter on the nights of April 7, April 28, 
and May 6 2008 UT with the 3.5-meter telescope at 
Apache Point Observatory (APO). We used a backside- 
illuminated SITe 2048x2048 CCD with 2x2 binning, re- 
sulting in a plate scale of 0.28"/pixel, and sub-framed 
to a field of view of 4.8' by 0.56' to decrease readout 
time. We applied typical overscan, bias, and fiat-field 
calibrations. For photometric reduction we used the 
standard IRAF task PHOT, with the aperture selected 
as a constant multiple of the Gaussian-fitted FWHM 
of each image to account for any variable seeing. We 
performed differential photometry with respect to the 
star USNO 1167-0208653 (2MASS ID 175252970) located 
at Ilh42ml2.08s, +26°46'07.45" J2000. This star has 
V=10.82 and color V-I=1.48, compared to Gl 436 which 
has V=10.68, and color V-I=1.70. In the error bar com- 
putation, we account for both standard noise from the 
photometry, as well as due to scintillation following equa- 
tion 10 of Dravins et al. (1998). Having obtained at least 
30 minutes of data on each side of the transit, we sub- 
tracted a linear fit for all data outside of transit vs. air- 
mass to account for any differential reddening. Resulting 
individual data points have errors ranging from 1.5 to 2.8 
mmag, which agrees with the rms of the residuals from 
the model fits, and a typical cadence of about 17 seconds. 
We have searched for correlated noise on th e timescale 
of ing ress and egress, via the technique of I Pont et all 
(|2006r ). but only find a statistically significant amount 
for the night of April 7, measured to be 0.11 mmag. The 
three transits are shown in Figure 1. 

We also carried out accompanying observations with 
the New Mexico State University (NMSU) 1-meter tele- 



scope at APO, in the V fiher on the night of April 7 2008 
UT, and in the I filter on the night of April 28 2008 UT. 
A 2048x2048 E2V CCD was used with 1x1 binning and 
sub-framing, resulting in a field of view of 8.0' square and 
a plate scale of 0.47"/pixel, and we applied the aforemen- 
tioned standard calibration and photometric extraction 
techniques. We performed ensemble photometry with 
respect to the USNO star that was used as the 3.5m ref- 
erence, as weh as BD-h27 2046 (V=10.64, V-I=0.44), and 
another star at Ilh42m00s, -|-26°45'56" J2000 (V=12.81, 
V-I=1.46). Resulting typical errors on individual points 
range from 3 to 5 mmag with a typical cadence of about 
12 seconds. 

The NMSU 1-meter telescope can also function as a 
robotic telescope, and is used intermittently to pho- 
tometrically monitor stars with known radial-velocity 
discovered planets to search for transits. A search of 
the 1-meter archives revealed that it observed Gl 436 
on the night of January 11 2005 UT, during which a 
transit should have occurred, according to the precise 
ephemeris for Gl 436b that is now available by incor- 
porating the many observed transits in 2007 and 2008. 
At the time, this 1-meter program depended on visual 
inspection of automatically generated photometry and 
plots. For this night, the plot had large temporal and 
brightness ranges, and thus the tiny transit was easily 
missed. However, now carefully inspecting the region 
constrained by the ephemeris, as well as re-performing 
the photometry to maximize signal-to-noise, we find a 
transit signature within a minute of that predicted by the 
ephemeris with reasonable width and depth, as shown in 
Figure 1. Individual data points have an error of about 
4 mmag, a cadence of 30 seconds, and we do not detect 
any correlated noise with any level of significance. 

We also conducted observations on the nights of April 
28 and May 13 2008 UT using a 24" telescope located 
at the Sommers-Bosch Observatory (SBO) on the Uni- 
versity of Colorado at Boulder campus, using an I filter. 
These observations also used a windowed chip and an 
exposure time to maximize signal-to-noise without satu- 
rating, and have comparable temporal resolution to the 
3.5m and Im telescopes due to a shorter readout time. As 
well, we used an unfiltered 11" telescope at Cloudcroft, 
NM (CC) with a SBIG ST-7E CCD and 2x2 binning on 
May 6 2008 UT, with a resulting cadence of about 25 
seconds. We have also gathered all the amateur data 
currently available on the system as compiled by Bruce 
Gary (http: / /brucegary net/ AXA/G J436 /gj436.htm) . 

3. MODELING AND DERIVATION OF PARAMETERS 

We use the JKTEBOP code (jSouthworth et all 
I2004al|bh to model all the tr ansit light curves m a con- 
sistent and uniform manner. ISouthworthI (|2008[ ) has re- 
cently performed an exhaustive analysis of fourteen tran- 
siting planets using the JKTEBOP code, and shows it 
compares well with results reported elsewhere. JKTE- 
BOP offers the advantage of incorporating a Levenberg- 
Marquardt optimization algorithm, improved limb dark- 
ening treatments, and extensive error analysis routines, 
which are critical for confirming any trends in the system. 

For each transit curve, we solved for the ratio of radii 
(k = Rp/Rj), the orbital inclination (i), the time of mid- 
transit (To), and a scale factor that defines the normal- 
ized value of the out-of-transit flux in the light curves. 



3 



In order to obtain reasonable results for the scale of the 
system for all data sets, t he sum of the radii (Rs + Rjp) 
was set to that found by iTorrei (120011 ). We also fixed 
the eccentricity to a value of 0.15 and the l o ngitu de of 
periastron to 34 3° as given by Demi ng et al.l (2007) and 
iMardlin^ (j2008f ). We used a q uadratic limb- darkening 
law with coefficients taken from IClaretl ()2000[ ) for Te/ / 
= 3500K, log(g) = 4.5, Vt = 2.0 km s-\ and [M/H] 
= 0.0, for the appropriate filters. In the case of the 
Sp itzer 8/xm data, we u sed the coefficients as determined 
bv lGillon etlHI (|2007a[) . From each fit, still assuming a 
constant sum of radii, we were thus also able to calculate 
the individual star and planet radii, as well as the depth 
and width of transit. In order to rule put any potential 
correlations in derived planet size and inclination, we 
then re-modeled all data with the same procedure, but 
also fixin g k, and thus the star and planet sizes, to that 
found bv iTorrei ()2007[ ). This generally produced simi- 
lar results, but for the noisier data sets achieved more 
consistent results. Parameters from both techniques are 
shown in Table 1. 

In order to obtain robust errors, we ran 1,000 Monte- 
Carlo simulations for each d ata set and perfor med a 
residual-permutation analysis ([Jenkins et al.l[200^ to in- 
vestigate temporally correlated noise. In both cases, 
the previously fixed parameters, as well as the limb- 
darkening coefficients, were allowed to vary so that their 
individual uncertainties would be taken into account in 
the derived parameter uncertainties. For each Monte 
Carlo simulation, random Gaussian noise with ampli- 
tude equal to the given error bars, or in the absence 
thereof the standard deviation of the residual scatter 
from the best-fit solution, was added to each data point 
and the curve re-fitted with random perturbations ap- 
plied to the initial parameter values. This ensured a de- 
tailed exploration of the parameter space and parameter 
correlations. However, this Monte Carlo technique will 
underestimate errors for certain parameters in the pres- 
ence of temporally correlated noise, which can result from 
trends in seeing, extinction, focus, or other atmospheric 
or telescope related phenomena (Southworth 2008). The 
residual-permutation method takes the residuals of the 
best-fit model, shifts them to the next data point, and 
finds a new solution. The residuals are shifted again, a 
new fit is found, and the process repeats as many times 
as there are datapoints. Thus, there is a distribution of 
fitted values similar to the Monte Carlo technique, but 
any temporal trends will have been propagated around 
the light curve, and thus taken into account. For our 
final errors we adopt the larger value found between the 
two methods, although for the majority of parameters 
and data sets the two methods agree quite well. 

In total we modeled 28 light curves, (16 professional 
and 12 amateur), covering 19 separate transit events over 
a baseline of nearly 3.3 years. 

4. TRANSIT TIMING AND ECLIPSE VARIATIONS 

Using the derived time of minima in Table 1 for all 
the data when allowing k to vary, we derive a new linear, 
error-weighted ephemeris of Tc(HJD) = 2454222.6164(1) 
-|- 2.643897(2)-E, where the parentheses indicate the 
amount of uncertainty in the last digit, and E is the epoch 
with E = the initial transit discovery of Gillon et al.l 
(|2007b[) . Using this ephemeris, we then compute an ob- 
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Fig. 2. — O-C diagram for all professional times of minima. 



served minus calculated (O-C) diagram for the time of 
transit center, as shown in Figure 2. We have currently 
excluded the amateur data from the plot due to much 
larger error bars, so that the high-precision data points 
can be seen clearly. We have examined the TTVs and 
various subsets thereof using a phase dispersion mini- 
mization technique ( StcUingwcrl 1978), but do not find 
any periods with statistical significance. Examining the 
best data, specifically the previously published data and 
our 3.5-meter observations, there is a standard deviation 
of 52 seconds. Assuming a sinusoidal TTV trend, we can 
then rule out any TTVs with amplitude greater than ~ 
1 minute. 

We have searched for any trends in derived inclina- 
tion, width, and depth of transit over time via error- 
weighted least-squares linear regression. In addition, we 
have also performed 10,000 Monte Carlo simulations for 
each fit, where gaussian noise with amplitude equal to 
each point's error bars was added in each iteration and 
the data re-fitted, with resulting la parameter distribu- 
tions giving robust errors. The two methods agree to 
within 1% for all values. As mentioned in §3, we mod- 
eled all the light curves by both allowing the ratio of radii 
to vary as well as fixing it, and thus we list the values for 
each set. Performing fits to all the data, we have a ten- 
tative detection of increasing inclination, transit width, 
and transit depth with time, as shown in Table 2. We 
present these fits with the actual data derived when fix- 
ing the radii in Figure 3. As a precaution against any 
bias being introduced by the much larger number of data 
points at later epochs, we decided to separately bin the 
2005, 2007, and 2008 data using an error-weighted mean, 
and re-fit the three resulting data points for each model- 
ing method. As shown in Table 2, the values agree very 
well with those derived when not binning the data. 

The trends are moderately dependent on the single 
2005 transit data point, which greatly extends the tem- 
poral baseline, and as such we are cautious about any 
claims. Resulting temporal trends when removing the 
2005 data point are also shown in Table 2. Although 
while removing the 2005 data point significantly weak- 
ens the claim of a variation of inclination with time, the 
trend of increasing width still holds. Also of interest is 
that at a rate of 0.120 deg yr^^, as derived from our 
fit to all the data fitted with a variable radius, the JK- 
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TABLE 2 

Trends in derived inclination, width, and depth of 
transit over time 



Radii 


Data Set 


deg yr ^ 


min yr ^ 


mmag yr ^ 


Variable 


All 
Binned 
No 2005 


0.120±0.062 
0.126±0.061 
0.092±0.099 


3.43±1.01 
3.53±0.97 
3.10±1.10 


0.28±0.16 
0.26±0.14 
0.29±0.17 


Fixed 


All 
Binned 
No 2005 


0.069±0.051 
0.071±0.050 
0.020±0.099 


2.36±0.84 
2.37±0.81 
1.68±1.29 


0.32±0.20 
0.32±0.19 
-0.01±0.42 
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Fig. 3. — Measured inclination, width, and depth of transit over 
time for all data, with the star and planet radii fixed. 

TEBOP program yields an increase in transit width of 
4.36 min yr~^, and depth of 0.544 mmag yr^^, which 
are in agreement with our observed trends, and thus are 
self-consistent. As well, the measured rate of inclination 
change is compatible with the ~0.1 deg yr~^ requ i red to 
make congruent the non- dete ction oflButler et al] (|2004f ) 
and the observed transit of iGillon et al.l (|2007b| ). Ex- 
tending the measurement baseline a couple years into 
the future will confirm or negate this result. 

5. DISCUSSION AND CONCLUSION 

We have presented a total of ten new transit light 
curves of Gl 436b, three of which come from the 3.5- 
meter telescope at APO, and one of which is from the 
NMSU 1-meter in January 2005. We have collected and 
uniformly modeled all available professional and amateur 
light curves, and searched for any trends in transit tim- 
ing, width of transit, and depth of transit variations. We 
find statistically significant, self-consistent trends that 



are compatible with the perturbation of Gl 436b by a 
planet with mass < 12 in a non-resonant orbit with 
semi-major axis < 0.08 AU. This conclu s ion is b ased on 
the numerical simulations of lRibas et al.l ()2008al . see Fig. 
1) who constrain the mass and semi-major axis of the 
theoretical second planet by examining which configura- 
tions could produce the observed orbital perturbations 
while still remaining undetected by the existing radial- 
velocity data. From our analysis, we infer a non-resonant 
orbit based on a lack of detected TTVs with amplitude 
> 1 minute. 

We stress that our measured trends are moderately 
dependent on our 2005 data, and thus subsequent high- 
precision observations over the next few years need to 
be carried out to confirm or refute this trend. If con- 
firmed, it would be strong evidence for the first extraso- 
lar planet discovered via orbital perturbations to a tran- 
siting planet. Also, we would like to note that although 
lAlonso et al.l (|2008f ) had previously limited the rate of in- 
clination change to 0.03±0.05 deg/yr, they did so only by 
measuring the change in width between the 2007 Spitzer 
observations and their own 2008 H-band data, which they 
found to be 0.5±1.2 minutes. Via Table 1, we find the 
difference in transit width between the two observations 
to be 1.5±1.4 minutes, which is in agreement with our 
derived inclination and width values, and is a more reli- 
able result due to using full model fits with proper limb- 
darkening coefficients. With respect to the the amateur 
observations, although they are numerous, the very small 
depth of the transit makes it a challenge for most small 
aperture systems, resulting in very large uncertainties in 
i and Tq. Also, while amateur observers are aware of the 
importance of precision timing, we of course cannot ex- 
amine each of their observing set-ups, and thus one must 
be aware of the possibility, although small, of systemic 
time offsets on a given night when interpreting their data. 
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TABLE 1 

Parameters derived from using the JKTEBOP code with Icr errors 



Epoch 


Source 


Filter'^ 


Tmin 


Inclination 


R* 


Rp 


Depth 


■Width 








(HJD-2450000) 


(deg) 


(R©) 


(Re) 


(mmag) 


(minutes) 




Ratio of Radii k Allowed to Vary 


-318 


NMbU Im 


■V 


3381. 85584±0. 00179 


86.02±0.23 


0.446±0.046 


6.19±4.81 


7.05±1.35 


47.0it7.1 





Gillon et al. (■2007b)'' 


V 


4222. 61617±0. 00060 


86.38±0.18 


0.463±0.016 


4.32±0.24 


6.98±0.43 


60.1±1.6 


1 


Shporer et al. ("20081'' 
ShDorer et al. ('2008')'' 
ShDorer et al. ('2008')'' 


None 


4225.26052±0.00089 


86.43±0.17 


0.463±0.015 


4.33±0.27 


7.16±0.81 


61.4±2.5 


1 


V 


4225.26050±0.00072 


86.35±0.17 


0.462±0.016 


4.47±0.25 


7.31±0.46 


59.2±1.9 


9 


R 


4246.41012±0.00079 


86.27±0.18 


0.456±0.014 


5.10±0.66 


9.07±0.87 


56.7±2.6 


22 


Gillon et al. f2007a') 


8^m 


4280.78219±0.00011 


86.34±0.16 


0.464±0.016 


4.23±0.16 


7.46±0.10 


59.8±1.1 


110 


Gregor Srdoc'' 


R 


4513.43393±0.00174 


86.10±0.24 


0.457±0.030 


5.04±2.91 


7.11±1.29 


49.8±6.5 


110 


Tonny 'Vanmunster^ 


R 


4513.44404±0.00247 


87.13±0.30 


0.461±0.016 


4.57±0.45 


9.63±1.67 


77.8±4.9 


112 


Bruce Gary^ 


R 


4518.72999±0.00278 


86.04±0.34 


0.443±0.087 


6.55±9.38 


8.78±2.37 


47.0±11.7 


113 


Gregor Srdoc'' 


R 


4521.37338±0.00130 


87.27±0.30 


0.459±0.015 


4.81±0.30 


11.06±1.13 


80.0±4.2 


115 


James Roe^ 


V 


4526.65995±0.00124 


86.03±0.27 


0.447±0.036 


6.05±3.34 


7.47±1.55 


48.4±9.0 


115 


Joao Gregorio" 


V 


4526.65972±0.00130 


87.09±0.28 


0.468±0.016 


3.80±0.25 


6.50±0.66 


77.1±4.0 


117 


Richard Schwartz^ 


V 


4531. 94399±0. 00222 


86.18±0.20 


0.461±0.015 


4.51±0.65 


6.37±1.67 


53.4±6.5 


118 


Alonso et al. (2008') 


H 


4534.59611±0.00014 


86.39±0.17 


0.463±0.016 


4.32±0.17 


7.74±0.11 


61.1±0.9 


127 


Manuel Mendez'' 


R 


4558.38849±0.00173 


86.60±0.24 


0.466±0.016 


3.99±0.33 


6.60±0.87 


66.6±5.6 


129 


NMSU Im 


V 


4563.67937±0.00257 


86.45±0.33 


0.467±0.019 


3.87±0.61 


6.34±1.12 


61.5±8.0 


129 


APO 3.5m 


V 


4563.67968±0.00051 


86.44±0.17 


0.459±0.015 


4.73±0.28 


8.60±0.44 


61.8±2.2 


132 


James Roe'' 


B 


4571.61844±0.00107 


88.60±0.62 


0.455±0.015 


5.24±0.24 


14.84±0.83 


95.5±3.6 


137 


NMSU Im 


I 


4584.83301±0.00117 


86.55±0.19 


0.449±0.015 


5.89±0.51 


14.61±1.72 


65.2±3.6 


137 


APO 3.5m 


V 


4584.83084±0.00035 


86.32±0.16 


0.464±0.015 


4.20±0.20 


6.36±0.25 


58.3±1.2 


137 


SBO 24" 


I 


4584.82868±0.00166 


86.63±0.21 


0.448±0.015 


5.95±0.50 


15.23±2.13 


67.3±3.8 


137 


Bruce Gary^ 


R 


4584.82876±0.00087 


86.51±0.18 


0.463±0.015 


4.32±0.25 


7.41±0.58 


64.0±2.4 


138 


Manuel Mendez^ 


R 


4587.47754±0.00170 


86.91±0.28 


0.462±0.016 


4.45±0.35 


8.83±1.13 


73.9±4.9 


140 


CC 11" 


None 


4592.76123±0.00140 


86.25±0.17 


0.463±0.015 


4.38±0.48 


6.64±1.03 


56.0±3.4 


140 


APO 3.5m 


V 


4592.76281±0.00084 


86.50±0.17 


0.465±0.015 


4.12±0.21 


6.71±0.36 


63.5±3.4 


140 


SBO 24" 


I 


4592.76202±0.00177 


86.55±0.26 


0.453±0.017 


5.43±0.84 


12.25±1.50 


65.3±4.8 


143 


SBO 24" 


I 


4600.69795±0.00118 


85.88±0.24 


0.425±0.067 


8.52±6.54 


6.75±1.08 


42.0±8.5 


146 


James Roe^ 


V 


4608.62470±0.00107 


86.32±0.23 


0.454±0.015 


5.31±0.67 


9.86±0.55 


58.4±6.3 




3.5m Data Combined 


V 




86.39±0.16 


0.463±0.015 


4.39±0.22 


7.25±0.31 


60.6±1.3 



Star and Planet Radii Fixed by Fixing k 



-318 


NMSU Im 


V 





Gillon et al. f2007b')'' 


V 


1 


ShDorer et al. ('2008')'' 


None 


1 


Shporer et al. ('2008')'' 
Shporer et al. ('2008')'' 


V 


9 


R 


22 


Gillon et al. f 2007a') 


8^m 


110 


Gregor Srdoc'' 


R 


110 


Tonny Vanmunster^ 


R 


112 


Bruce Gary'' 


R 


113 


Gregor Srdoc'' 


R 


115 


James Roe'' 


V 


115 


Joao Gregorio" 


V 


117 


Richard Schwartz" 


V 


118 


Alonso et al. (20081 


H 


127 


Manuel Mendez" 


R 


129 


NMSU Im 


V 


129 


APO 3.5m 


V 


132 


James Roe" 


B 


137 


NMSU Im 


I 


137 


APO 3.5m 


V 


137 


SBO 24" 


I 


137 


Bruce Gary" 


R 


138 


Manuel Mendez" 


R 


140 


CC 11" 


None 


140 


APO 3.5m 


V 


140 


SBO 24" 


I 


143 


SBO 24" 


I 


146 


James Roe" 


V 




3.5m Data Combined 


V 



3381.85596±0.00212 
4222.61617±0.00062 
4225.26049±0.00094 
4225.26050±0.00076 
4246.41009±0.00103 
4280.78219±0.00011 
4513.43416±0.00191 
4513.44424±0.00386 
4518.73038±0.00358 
4521.37312±0.00244 
4526.66055±0.00302 
4526.65996±0.00101 
4531.94392±0.00198 
4534.59610±0.00014 
4558.38809±0.00164 
4563.67966±0.00252 
4563.67971±0.00116 
4571.61831±0.00467 
4584.83373±0.00379 
4584.83084±0.00036 
4584.82787±0.00912 
4584.82874±0.00100 
4587.47761±0.00204 
4592.76119±0.00142 
4592.76248±0.00093 
4592.76090±0.00430 
4600.69668±0.00171 
4608.62542±0.00508 



86.15±0.17 
86.40±0.16 
86.45±0.19 
86.39±0.16 
86.37±0.17 
86.34±0.16 
86.19±0.18 
87.29±0.54 
86.20±0.25 
87.34±0.39 
86.14±0.27 
86.65±0.19 
86.26±0.26 
86.40±0.16 
86.51±0.21 
86.38±0.23 
86.53±0.18 
88.62±1.02 
86.84±0.78 
86.32±0.16 
86.64±0.73 
86.53±0.18 
87.00±0.31 
86.27±0.15 
86.47±0.17 
86.71±0.32 
86.08±0.17 
86.40±0.68 
86.43±0.16 



0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.017 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.017 
0.464±0.016 
0.464±0.015 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 
0.464±0.016 



4.23±0.28 
4.23±0.28 
4.23±0.30 
4.23±0.30 
4.23±0.29 
4.23±0.30 
4.23±0.28 
4.23±0.28 
4.23±0.31 
4.23±0.31 
4.23±0.28 
4.23±0.29 
4.23±0.31 
4.23±0.28 
4.23±0.23 
4.23±0.28 
4.23±0.28 
4.23±0.28 
4.23±0.27 
4.23±0.29 
4.23±0.28 
4.23±0.28 
4.23±0.29 
4.23±0.29 
4.23±0.29 
4.23±0.28 
4.23±0.30 
4.23±0.28 
4.23±0.28 



5.61±0.63 
6.74±0.69 
6.86±0.76 
6.65±0.72 
6.66±0.70 
7.47±1.00 
5.92±0.74 
8.20±1.11 
6.00±1.03 
8.32±1.26 
5.56±1.30 
7.42±0.93 
5.95±1.08 
7.42±0.98 
7.00±0.76 
6.63±0.79 
7.06±0.78 
9.25±1.27 
7.70±1.19 
6.45±0.62 
7.05±1.45 
7.07±0.76 
7.92±0.99 
6.24±0.66 
6.94±0.68 
7.57±1.06 
5.54±0.71 
6.53±2.03 
6.82±0.67 



52.5±3.2 
60.5±2.6 
61.8±3.7 
59.9±2.1 
59.6±3.6 
59.5±1.0 
54.0±4.3 
77.5±9.0 
55.7±7.3 
80.1±6.3 
51.7±9.1 
67.3±3.4 
56.1±7.5 
61.1±1.0 
64.2±4.7 
60.1±5.3 
64.1±3.2 
93.3±6.5 
72.6±14.5 
58.0±1.8 
67.8±17.1 
64.3±2.8 
75.0±5.3 
56.5±2.9 
62.6±2.9 
69.5±7.8 
50.4±3.6 
63.8±19.6 
61.7±2.7 



Note. — All errors are Icr 

" Amateur Observer with data obtained from Bruce Gary, ^http : //brucegary ■ net /AXA/ G J436/ gj 436 . htmj ' Data were digitized from published 
plot*^ Johnson-Cousins System 



